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ABSTRACT 

Few chronobiologists may be aware of the regression-to-the-mean 
(RTM) statistical artifact, even though it may have far-reaching influences 
on chronobiological data. With the aid of simulated measurements of the 
circadian rhythm phase of body temperature and a completely bogus stimulus 
(unicorn milk), we explain what RTM is and provide examples relevant to 
chronobiology. We show how RTM may lead to erroneous conclusions re- 
garding individual differences in phase responses lo rhythm disturbances and 
how it may appear as though unicorn milk has phase-shifting effects and can 
successfully treat some circadian rhythm disorders. Guidelines are provided 
to ensure RTM effects are minimized in chronobiological investigations. 
(Chronobiology International 1X(6). 1041-1053, 2001) 

INTRODUCTION 

A fundamental aspect of chronobiological research is the examination of how 
certain stimuli affect the timing of biological rhythms in individual organisms. 
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Intuitively, one would expect the influence of a stimulus to depend on when that 
stimulus is introduced with respect to the initial phase of an individual biological 
rhythm. This hypothesis has been corroborated across a wide range of species, 
dependent variables, rhythm period lengths, and types of stimuli. For example, it 
is known that bright light delays the human circadian rhythm of body temperature 
when it is introduced before an individual's temperature curve reaches its mini- 
mum value. Conversely, bright light advances the temperature rhythm when intro- 
duced in the hours immediately after an individual's temperature minimum (I). 
Melatonin has opposite effects on rhythm phase with respect to the time it is 
administered relative to the time of the temperature minimum (2). 

Other research work has shown that the degree to which human circadian 
rhythms react to rhythm disturbances (e.g., transmeridian travel) depends on 
individual differences in initial circadian phase before the disturbance takes 
place (3). Similarly, experimental work involving simulated shifts of the light- 
dark cycle has shown how the degree of entrainment is related to individual 
phase before the shift occurred (4). Tn addition, bright light has been found to 
"normalize" the phase of circadian rhythms for individuals who have unusu- 
ally delayed or advanced rhythms before the light "treatment" is prescribed 
(5,6). 

All the above experiments have one thing in common: The change in phase 
of a circadian rhythm in response to a certain stimulus is examined in relation 
to the original circadian rhythm phase measured before the stimulus is intro- 
duced. The initial rhythm phase data are either those from individuals, some of 
whom will be toward the extremes (tails) of the population distribution, or the 
average phase from a sample(s) of individuals selected asymmetrically from the 
population of interest (e.g., samples of "late peakers" or "early peakers"). A 
statisticcil artifact that can influence such investigations of the change in a mea- 
sured variable over time is called regression to the mean (RTM). While notes of 
caution concerning RTM have been communicated in fields related to chronobi- 
ology, such as general medicine (7), blood pressure monitoring (8), general epi- 
demiology (9), and psychiatry (10), to date the artifact has not been formally 
addressed in a chronobiological journal. 

Our first aim is to outline the characteristics of the RTM statistical arti- 
fact. Second, we discuss how the results of some typical chronobiological ex- 
periments may be influenced by RTM and so should be interpreted with ex- 
treme caution. Third, we outline the various coping strategies that can be 
employed to minimize the influence of RTM. Specifically, we show that sev- 
eral of the apparent phase-shifting characteristics of stimuli measured in rela- 
tion to initial phase can be found without any stimuli actually being introduced. 
As a tribute to the work of Cole (11), who first showed how circadian cycles 
can be erroneously detected in a "unicorn" (i.e., in completely random data), 
we give our "phantom" stimulus the label "unicorn milk." This substance rep- 
resents any stimulus that might be found to affect circadian rhythm phase, but 
is in reality an erroneous artifact of RTM. 
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DEFINITION OF REGRESSION TO THE MEAN 

The RTM effect has been described as a "chance finding masquerading as 
a real one" (7, p. 241). In essence, when a sample of subjects is investigated in 
a test-retest fashion, it will be extremely likely that subjects who score higher 
than average values on the initial test will show lower values, closer to the sam- 
ple mean, on a retest (10). Conversely, subjects who score particularly low values 
on the initial test will be likely to show an increase on the second test, and this 
increase is purely a statistical artifact. Regression to the mean is a mathematical 
law that "occurs with any variable that fluctuates within an individual, either 
genuinely or due to measurement error" (7, p. 241, italics added). Regression to 
the mean can occur irrespective of whether, prior to the second measurement, 
the subject is aware of the initial value in the pretest, 

One can demonstrate the RTM effect by simulating some data according to 
the methods of Yudkin and Stratton (7). One benefit of such a simulation is that 
it proves that RTM is a statistical law, independent from any human bias in data 
collection and interpretation. We have done such a simulation for some data 
relevant to chronobiology (Table I). First, we simulated typical phase markers 
of circadian rhythms for 45 phantom individuals. We chose the time of the tem- 
perature minimum T win as our phase marker and then simulated data that were 
normally distributed with a mean T min of 05:00 and a standard deviation of 0.9h. 



Table /. Details of Simulated Data 



A: 


B: "True" 






H: 1st 


F: 2nd 


Subject 


T 

1 www 


C: Error 


D: Error 


Measurement 


Measurement 


I 


5.53 


-2.03 


-1.25 


3.50 


4.28 


2 


5.10 


-0.44 


-0.02 


4.66 


5.08 


3 


5.23 


-1 .21 


1.10 


4.02 


6.33 


4 


4.78 


-0.94 


0.86 


3.84 


5.64 


5 


4.80 


0.89 


-0.96 


5.69 


3.84 


6 


6.12 


J. 85 


-1.03 


7.97 


5.09 


7 


6.24 


0.81 


-0.38 


7.05 


5.86 • 


8 


3.24 


-0.33 


0.20 


2 91 


3.44 


9 


5.47 


-1.73 


0.22 


3.75 


5.69 


10 


4.57 


-0.41 


2.09 


4.17 


6.66 




/ 


/ 


/ 


/ 


/ 


45 


4.56 


0.99 


-1.66 


5.56 


2.90 


Mean 


5.10 


-0.02 


0.12 


5.08 


5.23 


SD 


0.89 


1.05 


0.90 


1.45 


1.23 



"True" values ot 7* (li() (column B) were generated for 45 ' phantom" subjects 
(column A). Two estimates (columns E and F) of each subject's true T mm were 
made by adding random error components (columns C and D) to each subject's 
true value. 
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These parameters were chosen according to the population data collected by Baehr 
et al. (12) from normal healthy individuals. These data, shown in column B of Table 
1, represent the hypothetical "true values" of each individual's 7 min (Le„ the "real" 
values that hold no measurement error). 

We stress here that the true value of any variable for an individual can never 
be known with absolute certainty, but is merely estimated from measurements that 
are made in reality. Nevertheless, the more measurements of circadian rhythm 
phases that are obtained (i.e., the more 24h periods that are examined in each treat- 
ment condition of the experiment to form an average phase for each person in each 
experimental condition), the less measurement error there would be consequently. 
We stress that the reader should not think that the basic chronobiological tactic of 
taking multiple measurements over time within a biological rhythm cycle controls 
for RTM effects. Measurements are usually taken over time by chronobiologists for 
the ultimate calculation of biological rhythm parameters (e.g., phase); these mea- 
surements will show measurement error over repeated cycles. It is the test-retest 
error in rhythm parameter estimation that is our concern, not error in individual 
measurements made within a cycle. Therefore, our discussion of measurement error 
and RTM pertains most to research in which only one estimate of phase (one cycle 
of a biological rhythm) is obtained before and after any stimulus is introduced. 

To return to the data simulation, we then imitated the scenario of a re- 
searcher making two estimates of the true phase markers before and after expos- 
ing the individuals to our phantom chronobiologically active stimulus, unicorn 
milk (i.e., no effect at all of the proposed chronobiotic). Two new columns of 
data were formulated to represent the random error component that would be 
added to each subject's true value. Both these columns of error data were simu- 
lated with the same a priori stipulated characteristics: a mean of zero, a standard 
deviation of about lh, and distribution in a normal (Gaussian) fashion. The "er- 
ror" data in column C were then added to the true values to represent the initial 
phase estimate (column E). The error data in column D (with the same character- 
istics as in column C) were also added to the true values to give the poststimulus 
phase estimates (column F). Such a scenario is completely consistent with what 
happens in a real experiment; that is, the researcher can never measure the true 
value of anything. The researcher can merely estimate it since there is biological 
and instrument error associated with all quantitative measurements (13). 

The data from columns E and F, which represent pre- and poststimulus 7; )in 
values, are presented in Fig. I. One can see that the "spread" of T mn values is 
slightly less after the unicorn milk is given compared to the initial r n>in estimates. 
Individuals who had initial 7 mjn data that were unusually early or late are ex- 
pected to "move" closer to the average in the poststimulus test. This effect is 
completely due to RTM and is most apparent when the average of several i indi- 
viduals at the extremes of the distribution are calculated, as is the case in the 
first two examples in the next section. It should be noted that the standard devia- 
tion of the data from the whole sample does not have to differ between test and 
retest for RTM influences to be apparent. 
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Test 

Figure 1. A demonstration of regression to the mean with simulated data. One can see that the 
individuals with measurements furthest from the sample mean on the initial measurement seem to 
"regress" toward the mean on the second measurement, 

EXAMPLES OF REGRESSION TO THE MEAN 
IN CHRONOBIOLOGICAL RESEARCH 

Regression to the mean may have an impact on any study design that in- 
volves the examination of change in specific subsamples of a population or in 
relation to the range of initial values in the population. Regression to the mean 
will be found to some degree with any quantitative measurement (14-16). In- 
deed, several researchers from a wide range of backgrounds have been "caught 
out" by RTM. For examples relevant to medicine and psychiatry, see Refs. 17 
and 10, respectively. Regression to the mean may even explain why some amount 
of publication bias is inevitable in scientific journals since there is always judg- 
ment error between reviewers (18). We now provide three examples of studies, 
some real and some hypothetical, that are relevant to chronobiology and that 
may have been compromised by RTM influences. 

Treatment of Orcadian Rhythm Phase " Disorders" 
with Unicorn Milk 

There are individuals with marker phases of sleep and other circadian 
rhythms that are unusually early or late enough to cause social and medical prob- 
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lems. Lewy et al. (5) identified individuals as initially suffering from "advanced 
sleep phase syndrome." These people have particularly early times of retiring 
and waking and circadian rhythm acrophases relative to the general population. 
Lewy et al. (5) treated these subjects with bright light and found that sleep and 
circadian rhythm phases became later, closer to the "norm" for the general popu- 
lation. In a subsequent study, Lewy et al. (6) identified individuals with particu- 
larly late sleep and circadian rhythm phases (people suffering from "delayed 
sleep syndrome"). Subsequent treatment with light caused the rhythms to be 
phased at earlier times, bringing the individuals, again, closer to normal. While 
we are not suggesting that the results of Lewy et al. (5,6) were fully compro- 
mised by RTM, we can show that RTM may have confounding effects on any 
study that involves a treatment administered between two observations of some 
chronobiological variable taken from a single subject or sample. 

The data shown in Fig. 1 can also be used to simulate a study designed to 
examine whether a certain stimulus (in our case unicorn milk) can help "cure" 
unusually delayed or advanced rhythms, We ranked the T mm data measured in the 
pretest shown in Fig. 1 . We then identified the 5 people with the latest 7 min and 
the 5 people with the earliest 7 min in the sample. We then took their poststimulus 
data (also shown in Fig. I) and calculated the pre- and post- averages for the 
phased-delayed and phase-advanced people. These data are shown in Fig. 2. It 
can be seen that unicorn milk "caused" the average 7" min of the phase-delayed 
group to occur earlier and the average 7 min of the phase-advanced cohort to occur 
later. These changes were statistically significant (P < .001). 

Of course, nothing has happened between pre- and postmeasurements of 
phase. The apparent "normalizing" of phases is merely RTM of individuals who 
were selected from the "tails" of the population. An interesting aspect of RTM 




1 syndrome 

FYe-test Post-test 
Test 

Figure 2. A simulation of a study that examines whether a stimulus (unicorn milk) can bring the 
circadian rhythms ot "phase-delayed" and "phase-advanced" individuals back to "normal." The RTM 
effect makes if appear as though the unicorn milk is chronobiologically active; in fact, the individuals 
at the two exiremes of the population distribution are merely moving closer to the sample mean on 
the retest. 
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is that the more "unusual" the observations are (i.e., the more the observations 
can be placed in the extreme tails of the population distribution), the greater will 
be the RTM effect (7). This would mean that the more severe the individuals in 
a sample are suffering from either advanced or delayed sleep phase syndrome, 
the greater will be the apparent movement of the sample mean back to normal 
on a second measurement. Consequently, the greater will be the apparent effec- 
tiveness of the treatment. It should be noted thai the addition of "healthy con- 
trols" and observing no response to the treatment would not mitigate any RTM 
influences. These subjects would be closer to the population mean at the outset, 
so RTM influences would be minimal on the rctest. 



Individual Differences in the Responses to a Rhythm Disturbance 

In the article, "Phase Shift in Temperature Rhythm After Transmcridian 
Flight as Related to Pre-f light Phase Angle," Colquhoun (3) found that individu- 
als with a particularly late initial phase of the body temperature rhythm showed 
a greater shift in rhythm phase following an eastward flight across eight time 
zones than those individuals with a particularly early initial phase. It was stated 
that the preflight individual differences in phase were "abolished" by the flight 
since the greater shift by the late peakers was apparent throughout all the post- 
flight days. We can show that RTM can mediate findings similar to those of 
Colquhoun (3). 

Figure 3 shows the results of a simulation of Colquhoun's study (3). In our 
simulation, two more estimates of r min , in addition to the ones already made, 
were simulated and combined with the data shown in Fig. 1 to give four repeated 




Figure 3. A simulation of a study that examines the change in phase of a circadian rhythm 
following a flight across eight time zones in relation lo (he initial phase before the flight. Those 
individuals with a particularly late phase show a greater change in phase in response to the trans- 
meridian travel. This is not due to any inherent differences between individuals, but merely the 
influence of regression to the mean (RTM). 
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measurements of 7* niin for the 45 phantom subjects. The first measurement repre- 
sented the preflight estimates of phase to which all the postflight T min were com- 
pared (preflight 7",™ -zero in Fig. 3) and which were used to designate subjects 
as late peakers or early peakers* Simulated phase shifts of 6h, 7h, and 8h were 
added to the three respective postf light measurements of phase to simulate grad- 
ual adjustment of circadian rhythms over the 3 postflight days [Colquhoun (3) 
measured T min over 12 days postflight, but we simulated only 3 days for clarity). 
It is stressed that the enforced phase shifts were, in reality, equal for all our 
phantom subjects. The sample was divided into two groups of late peakers and 
early peakers according to the methods of Colquhoun (3). 

It can be seen in Fig. 3 that the late peakers shifted on average more than 
the early peakers, in agreement with the original Figs. 1 and 2 that were pre- 
sented in the article by Colquhoun (3). This difference between groups is still 
evident several days postflight, also in agreement with Colquhoun (3). It is 
stressed that the effect in our simulation is completely due to RTM, rather than 
any individual characteristic that can influence the degree of phase shift, because 
the same degree of phase shift (total of 8h) was simulated for all individuals that 
were "studied/' The apparent difference in phase shift between late peakers and 
early peakers does not get larger over the postflight days since it is known by 
statisticians that RTM effects are reduced with increased measurements (7). In 
other words, the RTM effect occurred mostly between the preflight and day 1 
postflight measurements and then did not influence the data, so the RTM effect 
remained throughout the subsequent postflight tests. 

Phase-Shifting Effects of Unicorn Milk 

We can also show that unicorn milk has phase-shifting effects; giving uni- 
corn milk at night a few hours before 7 min leads to phase advances of T min on 
average and giving unicorn milk in the early morning a few hours after 7 mir 
generally causes phase delays. 

We simulated an experiment in which unicorn milk was randomly adminis- 
tered to the 45 individuals at two times of day. Unicorn milk was given at either 
03:00, before the usual r mil1 , or at 07:00, after the conventional T mil] . Therefore, 

22 subjects obtained the unicorn milk at 03:00, before the mean initial r min , and 

23 subjects obtained the unicorn milk at 07:00, after the mean initial T nm . The 
time of administration of unicorn milk relative to each individual's 7 mm was plot- 
ted on the abscissa of Fig. 4. Phase changes were calculated by subtracting post- 
stimulus phase from prestimulus phase, and these data were plotted on the ordi- 
nate to form a "partial" phase-response curve for unicorn milk (Fig. 4). Phase 
delays are shown as negative values, and phase advances are shown as positive 
values, as is the convention for phase-response curves (19). 

One can see from Fig. 4 that there is some amount of random variation 
present (as in any experiment). Nevertheless, it is clear that there is a tendency 
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4 




-6 -4 -2 0 2 4 6 

Time of stimulus relative to Tmin (h) 



Figure 4. A simulation of a study that examines whether a stimulus (unicorn milk) can phase 
delay or phase advance circadian rhythms depending on when that stimulus is introduced in rela- 
tion to the initial phase. Because the initial phase is involved in the calculation of both ,r and v 
data, some relation between x and v would be expected purely due to the chance influences of 
regression to the mean (RTM). 

for phase changes before the initial 7 min to be positive (advances) and a tendency 
for phase changes after the initial T mln to be negative (delays). A common finding 
also shown in Fig. 4 is a "crossover point" between phase advances and delays 
centered at circadian time zero (the average T inm for the sample, 05:00 in our 
case). Although the phase-response curve as a description of the data is often the 
end point in chronobiological research, one can also formally analyze such 
changes in phase with respect to circadian time in various ways. For example, a 
trend line can be fitted to the data as in Fig. 4, and linear or nonlinear regression 
models can be employed (20). The Pearson's product moment correlation coeffi- 
cient for the data in Fig. 1 is -0.53 (P < .001). Alternatively, the phase changes 
before and after circadian time zero (T min ) can be compared by employing a 
hypothesis test, the choice of analysis of Campbell and Murphy (21), for exam- 
ple, A two-sample / test for the data in Fig. 1 suggests that the mean phase 
change of +0.4h before 7 min is significantly different from the mean phase change 
of -0.7h when unicorn milk is administered after 7 min (P < .02). Both these 
changes are significantly different from zero when analyzed with one-sample t 
tests (P< .05). The results of these hypothesis tests serve to illustrate the "gar- 
bage in, garbage out" mantra in statistics since all they confirm is that RTM 
influences are present in the data, rather than any chronobiotic systematic effect. 

How can random data with no stimulus-mediated difference between initial 
7 min and final 7 min appear to form a partial phase-response curve? The important 
point in Fig. 4 is that initial T inm is involved in the calculation of both the x-axis 
data (time of stimulus minus initial T min ) and the v-axis data (initial 7 min minus 
final 7 min ) on the phase-response curve. This means that the x and y data are 
already mathematically related, irrespective of the presence of any physiological 
relationship. Therefore, some degree of negative relationship (in the present ex- 
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ample) between a and y is guaranteed because of regression to the mean (22). For 
a fuller discussion of this topic and more examples related lo clinical research in 
particular, readers are referred to the work of Blomqvist (23). Regression-to-the- 
mean influences have been discussed by such authors as Blomqvist (23), even 
though the investigation of change relative to initial value may be a relatively 
small issue in a clinical study. The fact that, in chronobiological research, such 
an investigation is very likely to be the major issue serves to confirm the impor- 
tance of considering RTM influences. 

We stress that it cannot be claimed that every phase-response curve involv- 
ing separate subjects for each data point is .compromised by RTM since some 
phase-response curves show opposite trends from those expected due to the arti- 
fact (e.g., bright light). Moreover, some phase-response curves show a "dead 
zone" away from T mm , which is associated with zero phase change. The RTM 
could not explain such a characteristic of the phase-response curve. In this re- 
spect, it is stressed that the phase-response curve in Fig, 4 is merely a "partial" 
phase-response curve derived from only a few reference times of stimulus (03:00 
and 07:00) and did not assess whether a "dead zone" was present. Nevertheless, 
such plots can still be found in the literature (20). 



REGRESSION TO THE MEAN AND MEASUREMENT ERROR 

An important characteristic of the regression to the mean artifact is that the 
higher the correlation between the phase estimates before and after the stimulus, 
the less will be the influence of the artifact (15). This correlation between pre- 
and postestimates of phase is governed by the measurement repeatability of the 
particular variable. Thankfully, at least for body temperature, a moderate- to-high 
correlation (r = 0,78) between test-retest estimates of phase has been reported 
(12). Nevertheless, one can be absolutely certain that there would be no influ- 
ence of the artifact only if the repeatability correlation between phase estimates 
is unity (15). Since this scenario is almost impossible with physiological data 
measured on a quantitative scale (13), it is important that researchers ensure 
the artifact is controlled for in any chronobiological investigation. Streiner (10) 
proposed that researchers could use the following equation either to predict the 
magnitude of regression to the mean prior to the experiment taking place or to 
establish ex post facto how much a particular finding could be explained by 
regression to the mean: 

T 2 = X+r{T, -X) 

where 7, is the first measurement of a person or the mean of a sample, T 2 is the 
second measurement of a person or mean of a sample, X is the population mean, 
and r is the test-retest repeatability coefficient for the variable of interest Again, 
the above equation dictates that regression to the mean would be worst in scenar- 
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ios in which the subjects are extreme relative to the population mean and the 
test-retest correlation is relatively low. 

CONTROLLING FOR REGRESSION TO THE MEAN 
IN CHRONOJBIOLOGICAL RESEARCH 

To control for RTM, Campbell and Makin (16) and Bland and Altman (22) 
suggested that the average of prestitnulus and poststimulus data (initial and final 
T m \„), rather than initial T im(i per se, should be used as the baseline (prestimulus) 
data. In the case of Fig. 4, the *-axis data would be the reference time of unicorn 
milk minus the average of the two T mn values (pre- and poststimulus). While this 
control method for regression to the mean would remove the artifact from the 
analysis, the biological significance of a relationship between change in value 
and the average of final and initial values is difficult to interpret (24). Indeed, 
chronobiological research is unique in that the timing of the initial phase has 
intrinsic meaning. Taking the average of initial and final phases may compromise 
the whole interpretation of a phase-response curve. 

A preferable, but also more tedious, coping strategy for RTM is to take 
multiple measurements of the initial phase since RTM can be reduced by using 
the mean of several initial measurements (7). This tactic makes sense since it is 
known that the more repeated measurements are obtained, the closer will be the 
average of those measurements to the true value (13). Davis (9) showed that RTM 
is negligible when the average of four repeated measurements is used as the data 
involved in the analysis. Obviously, taking four measurements of phase prior to 
the introduction of any stimulus is inconvenient, but would be necessary to be 
confident that RTM influences are absent in certain research designs. 

Thankfully, the detrimental influence of RTM on the internal validity of a 
study can be controlled to some extent by the research design. Data drawn from 
double-blind experiments involving counterbalanced "repeated measures" would 
not be affected by the artifact since the same subjects have been examined re- 
peatedly under counterbalanced administration of placebo or stimulus. In this 
respect, the subjects are acting as their own controls. This research design is 
different from that examined in the above simulation of a phase-response curve 
since our stimulus was introduced only once to each individual at different times 
of day. We stress the importance of an adequate placebo-controlled treatment or 
group in any study that examines phase responses to a certain stimulus. 

CONCLUSIONS 

Any research that involves the measurement of pre- and postcircadian 
rhythm phase and the assessment of change in phase in relation to initial phase 
could be affected by RTM. The degree of influence on research conclusions 
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would depend on the measurement error that is present and the design of the 
study. Chronobiologists are asked to be wary of RTM effects in interpreting past 
research work and in designing their own studies. The equation cited in this 
article could be used to predict the magnitude of RTM influence in past or pro- 
posed single-group pre/posttest research designs. Ideally, the chronobiologist 
should take multiple measurements of initial phase to control for RTM influ- 
ences, especially if the repeatability of phase estimates is known to be relatively 
low. Inclusion of control groups or conditions administered in a counterbalanced 
order is encouraged in any research design. 
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